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Abstract 



We investigate some issues concerning the zero-momentum four-point renor- 
malized coupling constant g in the symmetric phase of O(N) models, and the 
corresponding Callan-Symanzik /3-function. 

In the framework of the 1/N expansion we show that the Callan-Symanzik 
/3-function is non-analytic at its zero, i.e. at the fixed-point value g* of g. 

This fact calls for a check of the actual accuracy of the determination 
of g* from the resummation of the d = 3 perturbative g-expansion, which is 
usually performed assuming the analyticity of the /3-function. Two alternative 
approaches are exploited. 

We extend the e-expansion of g* to 0(e 4 ). Quite accurate estimates of 
g* are obtained by an analysis that exploits the analytic behavior of g* as a 
function of d and the known values of g* for lower-dimensional O(iV) models, 
i.e. for d = 2,1,0. 

Accurate estimates of g* are also obtained by a reanalysis of the strong- 
coupling expansion of the lattice iV-vector model allowing for the leading 
confluent singularity. 

The agreement among the g-, e-, and strong-coupling expansion results 
is good for all values of N. However, at N = 0, 1, e- and strong-coupling 
expansion favor values of g* which are slightly lower than those obtained 
by the resummation of the g-expansion assuming analyticity in the Callan- 
Symanzik /3-function. 
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1/iV-expansion, e-expansion, Strong-coupling expansion. 
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I. INTRODUCTION 



The renormalization-group theory of critical phenomena provides a description of statis- 
tical models in the neighbourhood of the critical point. For O(N) models calculations are 
based on the 4 -field theory defined by the action 



S = / d d x 



(1) 



A strategy, which has been largely employed in the study of the symmetric phase, relies 
on a perturbative expansion in powers of the zero-momentum four-point renormalized cou- 
pling constant g performed at fixed dimension d = 3 Q. This perturbative expansion is 
asymptotic; nonetheless accurate results can be obtained by resummations exploiting its 
Borel summability and the knowledge of the large-order behavior. As general references 
on the g-expansion method see for instance Refs. This technique has led to accurate 

estimates of the critical exponents. 

An important quantity entering the calculation of universal quantities is the fixed-point 
value of g, i.e. the zero of the corresponding Callan-Symanzik /3-function. In the critical 
region, the bare coupling constant go becomes infinite on the scale fixed by the correlation 
length, whereas the zero-momentum four-point renormalized coupling approaches a finite 
non-zero limit g* at criticality. Accurate calculations of g* have been done by analyzing 
the perturbative expansion of the Callan-Symanzik /3-function [|5|-|T0H (known to 0(g 7 )), 
using the known results for its large-order behavior. The best determinations of g* have 
been apparently obtained by Le Guillou and Zinn- Justin M by making some additional 
assumptions on the analytic properties of the Borel transform. Such additional assumptions 
have been questioned by Nickel [|TT|,0], who argued the presence of confluent singularities 
in the /3-function at its zero, which may complicate the analytic structure of the Borel 
transform. 

This issue needs a non-perturbative analysis to be clarified, furthermore an analytic 
approach is required in order to understand the nature of the singularities. These features 
can be realized in the framework of the 1/N expansion. In this paper we analyze the Callan- 
Symanzik /5-function computed in Ref. ||12|| to Q(l/N) (i.e. the next-to- leading order) . While 
the leading order is analytic, the 0(1/N) term shows the presence of confluent singularities 
at the zero of the /3-f unction for all 2 < d < 4. Moreover a phenomenon analogous to the 
Abe-Hikami anomaly for the specific- heat [HJ emerges at the special dimensions d = 4 — 2/n 
(where n is an integer number), thus including the interesting case d = 3. 

In the analysis of Ref. pi confluent singularities at the zero of the /3-function may cause 
a slow convergence to the correct fixed-point value of g ||11|| . The apparent stability of the 



results when analyzing a finite number of terms of the perturbative expansion may then not 
provide a reliable indication of the uncertainty of the overall estimate. Confluent singularities 
represent a source of systematic error for the procedure used in Ref. §. A more general 
analysis explicitly allowing for the presence of confluent singularities would slightly change 
the value of g* for small values of N (although not excluding the values obtained in Ref. ||) 
and consequently the values of the critical exponents |J . It is therefore important to exploit 
other approaches to the study of O(iV) models, which can provide a check of the estimates 
of g* from the resummations of the perturbative ^-expansion, and of their actual accuracy. 
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An alternative field-theoretic strategy is the expansion in powers of e = 4 — d |14| . An 



important advantage of the e-expansion is the possibility of working directly at criticality. 
This allows us to go from one phase to another using the same framework. In order to get 
estimates at e = 1, the e-expansion requires eventually a resummation which is usually per- 
formed assuming its Borel summability. Relatively long series of the critical exponents have 
been calculated and their analysis has led to estimates which are in substantial agreement 
with those determined by the ^-expansion (see e.g. Ref. P~5fl). 



The fixed-point value of the zero-momentum four-point renormalized coupling is known 
only to 0(e 2 ) |TE], thus not allowing a real check of the value of g* obtained by the g- 
expansion. In this paper we extend this calculation to 0(e 4 ). The analysis of the 0(e 4 ) 
series of g* provides already good estimates with an apparent uncertainty of approximately 
6% for small values of N. A considerable improvement in the analysis of the e-expansion 
of g* is achieved using the known values of g* for lower-dimensional O(N) models, i.e. for 
d = 2,1,0. The key point is that g* is expected to be analytic and quite smooth in the 
domain < d < 4 (thus < e < 4). This can be indeed verified in the large- iV limit to 
0(1/N) using the results of Ref. ||12|| . Generalizing the technique presented in Ref. ||15|| , we 
perform a polynomial interpolation among the values of d where g* is known (d = 0, 1) or 
for which good estimates are available (d = 2, especially for N = 0, 1 by strong-coupling 
calculations), and then analyze the series of the difference. This procedure leads to much 
more accurate estimates of g*, which are consistent with those obtained by the direct analysis 
of the original e-series, but with an apparent error of approximately one per cent (~ 0.5% 
for N = 1). The agreement with the g-expansion estimates is good for all values of N. 
However it is worth anticipating that the results for N = 0, 1 turn out slightly lower than 
the estimates given by Le Guillou and Zinn- Justin || , thus favouring the results of the more 
general analysis done by Nickel ||. 

Another approach which has been widely used in the study of critical phenomena is 
based on lattice formulations of the theory. Two main techniques have been exploited in 
this context: high- and low-temperature expansions and Monte Carlo simulations. In the 
symmetric phase the fixed-point value of the zero-momentum four-point renormalized cou- 
pling can be estimated by analyzing the strong-coupling expansion of the quantities entering 



its definition |17-20,O,Elhp3]. Studies based on Monte Carlo simulations can be found in 



Refs. p4^pT|. The agreement with the field-theoretic estimates is substantially good, but 
small discrepancies have been observed, especially for small values of N |r2| , p2| , |30| , |3ll . We 
will show that such residual discrepancies disappear (or are largely reduced) when the lead- 
ing effects of the confluent singularities are properly taken into account in the analysis of 
the strong-coupling expansion [3"2,3"3",23j and in the analysis of the Monte Carlo data. 
The paper is organized as follows: 

In Sec. H we introduce some definitions and notations used in the paper. 

In Sec. [TTT]we analyze the Callan-Symanzik /5-function calculated to 0(1/N) in Ref. ||12|| , 
and discuss the presence of confluent singularities at its zero. 

In Sec. [TV] we analyze the expansion to 0(e 4 ) of the fixed-point value of the zero- 
momentum renormalized coupling constant. 

Sec. [V] is dedicated to a reanalysis of the strong-coupling expansion which allows for the 
confluent singularities present in the lattice approach. 

In Sec. [VT] some conclusions are drawn. 
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In App. |A| we give some technical details on the analysis of the analiticity properties of 
the Callan-Symanzik /3-function to 0(1/N). 

In App. JTJ we present the perturbative calculation to 0(e 4 ) of g*. 



II. DEFINITIONS AND NOTATIONS 



Let us introduce a few definitions and notations. If we set 

(r(0)^(x)) = S^G^(x), (2) 
{F{U)<P(x)<P{y)<j> 6 {z)) e = + S^S* + 5^5^)G^(x, y, z), (3) 

the zero-momentum four-point renormalized coupling g is defined as 

/ dxdydz G <y4 \x,y, z) 
9= £ d [fdx G( 2 )(x)] 2 ' 



where £ is the second-moment correlation length 

2 _ 1 Jdxx 2 G^(x) 
* 2d fdxGW(x) ' 1 J 

The normalization in Eq. ([|) is such that in perturbation theory g = go/m^ + 0(<7o/ m o e )- 

For convenience in the paper we will also introduce other definitions. First of all we 
consider the rescaled coupling 

1 + / d\ 

Unlike g* which is of order (4 — d) = e for d — > 4, the fixed-point value of g is 0(1), due 
to the factor multiplying g in its definition. Moreover g has the property that for N — > oo, 
g* — > 1 for any dimension d. For d = 3, it coincides with the coupling which is usually used 
in the analysis of the perturbative expansions in fixed dimension d = 3 |5||| . In the large- A 
limit another definition is also useful: 

9 = — ■ (7) 
Finally a fourth definition is quite common in the literature: 

/ S ^ (8, 
The quantity / is naturally defined in terms of 

X = J dx (0(0) (9) 
X 4 = / dxdycfe (0(0) • <t>{x)(j){y) ■ <f>(z)) c , (10) 
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as 



It is well-known that the model defined in Eq. ([!]) is equivalent at criticality to the 
iV-vector model defined in the continuum by 

S=^L [ d d x d„s{x) ■ d^s(x) (12) 



where s(x) ■ s(x) = 1. 

On the lattice one can consider any discretization with the (formal) continuum limit given 
by Eq. fll2T) . Here we will study the theory with nearest-neighbour interactions defined by 

S L = -f3NY,s(x)-s(y), (13) 

where the sum extends over all lattice links (xy). In this case the renormalized coupling 
constant is given by the previous formulae with the obvious substitution — > s. Of course, 
on the lattice, the integrals are replaced by sums over the lattice points. 



III. CRITICAL-POINT NON-ANALYTICITY IN THE LARGE- LIMIT 



An important controversial issue in the field-theory method at fixed dimension d = 3 
is the presence of non-analyticities at the critical point g*. The question was raised long 



ago by Nickel [ITJ who gave a simple argument to show that non-analytic terms should in 
principle be present in the /3-function. The same argument applies also to other series, like 
those defining the critical exponents: any quantity should be expected to be singular at the 
critical point. 

To understand the problem, let us consider the four-point renormalized coupling g as a 
function of the temperature T. For T — > T c we can write down an expansion of the form 



9 = 9" 



1 + ai (T - T c ) + a 2 (T - T c ) 2 + . . . + 6i(T - T C ) A + b 2 (T - T C ) 2A + . . . + 



ci(T - T C ) A+1 + . . . + di(T - T C ) A2 + . . . + ei(T 



T C ) A3 + . . . 



(14) 



where A, A 2 , ... are subleading exponents. We expect on general grounds that a\ = a 2 = 
— ... — 0. Indeed these analytic corrections arise from the non-linearity of the scaling 
fields and their effect can be eliminated in the Green's functions by an appropriate change 
of variables [34"]. For dimensionless renormalization-group invariant quantities such as g, the 



leading term is universal and therefore independent of the scaling fields, so that no analytic 
term can be generated. We will explicitly show their absence in the following computation. 
Notice that analytic correction factors to the singular correction terms are generally present, 
and therefore the constants q in Eq. ([14]) are expected to be nonzero. Moreover we mention 
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that, in general, the correction terms can include powers of the critical exponents^ 1 — a 
and 7 [0] (associated with backgrounds). 

Starting from Eq. ( |HD it is easy to compute the /3-function: 

^ = ~dM = dM/dTdT- (15) 

Since the mass gap M scales as 

M ~ (T — T c y [l + a l (T- T c ) + ... +b 1 (T - T C ) A + ...] , (16) 

we obtain the following expansion: 

(3(g) = ai (g* -g) + a 2 (g* - gf + . . . + fa(g* - g)* + (3 2 (g* -g)* + ...+ 

li(g*-g) 1+ * + ... + 6 1 (g*-g)% + ... + ( 1 (g*-g)% + .... (17) 

It is easy to verify the well-known fact that a\ = —A/v = —uo and that, if a\ — a 2 = . . . = 
in Eq. (0), then (3i = fa = ■ ■ ■ = 0. Eqs. (|i~4f) and dl7|) are the expressions which 
are expected on the basis of the standard renormalization-group picture. They express the 
asymptotic behaviour for generic models and indeed this is the supposed behaviour in lattice 
theories. However for the continuum A0 4 a much simpler expansion is often conjectured?] 
38| , |39|| . First of all one assumes that the continuum theory couples only to one subleading 



scaling field, the operator associated to the exponent A. As a consequence in the expansion 
(|14]) no corrections with exponents A 2 , A 3 , . . ., should appear. In Eq. (0) this conjecture 
implies for instance that 5i = Ci = 0, i.e. no terms with exponents Aj/A are present. The 
second claim is that T — T c = — m^ c is a scaling field in the Wilson renormalization- 
group sense fHfl . Therefore in Eq. (|TJ) no term with exponent hA + k, k > should appear. 



Correspondingly in Eq. (|T7|) terms of the form h + k/A would be absent. As a consequence 
(3(g) would be analytic. Of course, these two hypotheses would also prove that also other 
expansions in g, like the series for the critical exponents, would be analytic at the critical 
point. 

In order to understand the validity of these conjectures, one must compute the /3-function 
non-perturbatively. The only case in which this is possible is the large- iV limit, i.e. in the 
framework of the 1/N expansion. 



In the following we will analyze the /3-function computed in Ref. |L2] for N — >• oo and 
we will show explicitly that the conjecture mentioned above is incorrect: non-analytic terms 
are indeed present. However, at the order of 1/N we are working, we will be unable to 



1 In the following large- ./V analysis we will not consider these terms. Since for large values of N, 
7/V = (1 — a)/v = 2, their contributions mix with the analytic background, so that their inclusion 
would not change our main conclusions. 



2 For a critique of these conjectures from the point of view of the renormalization group a la Wilson, 



see Ref. [35|, Sec. 5.2 and App. E of Ref. 36]. See also Ref. [37] for a more recent discussion in the 
same framework. 
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distinguish which hypothesis is false. Indeed since, at N = oo, 1 + 1/A = A 2 /A = A 3 /A, 
we will only be able to verify that at least one of the corresponding terms is present but 
not to prove the presence of all of them. This problem can only be solved by computing the 
next two orders in 1/N. 

The starting point is the /^-function expanded in powers of 1/N: 



m = M *L = /,<»>© + I,«® + (i,) 



(18) 



where, for convenience, we have introduced the coupling g defined by Eq. (|7p. The two 
functions /3<®(g) and (3^{g) were computed in Ref. ||12|| . One has 



P 



(0), 



(d-A)g 1- 



(19) 



and 

P {1) {g) 



9 1 



(d - 3)2 d ' 1 /3 + - d (d - If (d - A + p g) 2 J 



d d u 



+2 



d d u 

(27T) a 



Pod + d - A 



Pod 



l + gU(u)Y 
d d u {p g + d-Af 
dJ (2nY (l+gU(u)) 2 



l + gll(u) 
1 



(2ir) d [1 + gU(u)} 2 (A + u 2 ) 
3 / 4 - 2 



+ 



d-1 



1 + U 2 ) 3 (1 + M 2 )(4 + M 2 ) \1 + u 2 A + u 2 , 



-4 
-4 



d d u 
(2tt) 



/3 ^ + rf-4 /5 ^+f-2 



;i + £Ih») d (l + ^n( M )) 2 J (4 + u 2 ) 



A P g(P g + d — A) 



/3 £ + d - 4 



(27r) d (l + ^n(w)) 2 (4 + w 2 ) 2 
Here g^ is the critical value of g for N = oo 

2(An) d / 2 



P g + d 
3 

1 +M 5 



1 + 



/?o is defined by 



and 



ft 



r(2-d/2)' 



4-d 



4e 



(2tt)V + 1 



(p + u) 2 + 1 p 2 + 1 



(1 + u 2 ) 3 (l + u 2 )(A + u 2 ) \l + u 2 A + u 2 



(20) 



(21) 



(22) 



(23) 



The leading term, Eq. (0), is clearly analytic. In App. [A] we study the behaviour of P^(g) 
for g -►</£,. Setting 



6 = 



(24) 
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for 2 < d < 4, we find 



analytic terms + A(d)Q i - d + . . . , (25) 



where 



•ly+i^d 2 d " 4 (M 4 - 112rf 3 + 428rf 2 -512rf+192) nN d 
U ~ ~ \c) 8d(Q-d)(A-d) sin(27r/(4-d))' ( ' 



and 



0F r(d/2-i) 

2 r((d-i)/2)' { ] 

2 

^ = (47r)^r(d/2)- (28) 

This expansion is not valid whenever d = 4 — 2/n, n e IN, as in this case A(d) diverges. For 
these values of the dimension one finds 

3^ (a) 

H w = analytic terms + B(n) 9 n log 9 + . . . (29) 

where 

In particular for d = 3, which corresponds to n = 2, we have 



B(2) = — (31) 
12tt 5 



Let us now interpret the results. Let us consider first the case of generic d / 4 - 2/n. In 



the large- N limit the smallest exponents (in the region 2 < d < 4) are j|I 



A = i-| + 0(l/JV), (32) 
A 2 = A 3 = ^2+0(1/JV). (33) 



At iV = oo A 2 and A 3 are degenerate independently of d |41| . This degeneracy should be 
lifted only at the next orderQ, i.e. 0(1/ N). To leading order in 1/N one has 

A 2 A 3 1 2 

X = X = 1 + A = 4^r < 34 > 



The corresponding scaling operators are linear combinations of 0\ = (toq + ^go<fi )(V0) and 
2 = (V 2 </>) 2 0. 
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Eq. (p5| ) shows that non- analytic terms are present. But because of the degeneracy 
they cannot be distinguished. Correspondingly in Eq. ([n|) one has 

7i + h + Ci ^ 0. (35) 

The evaluation to 0(1/N) of the constant 71, S± and (j requires an 0(1/N 3 ) computation. So 
we cannot check if some of them are zero. Notice that, as expected, no term with exponent 
1/A appears in Eq. (|25|) in agreement with the argument we presented at the beginning. 

Let us now consider d = 3 (an analogous argument applies to any special dimension d = 
4— 2/n). In this case the interpretation is more difficult and we have a phenomenon analogous 
to the Abe-Hikami anomaly for the specific heat [U| . The origin is an additional degeneracy 



of the exponents for iV — > 00: the non-analytic terms with exponents A 2 /A, A 3 /A and 
1 + 1/A become degenerate with the analytic term with exponent 2. This degeneracy 
causes the appearance of the logarithmic term in Eq. ( p9|) and has the consequence that the 
coefficients of the non-analytic terms are of 0(1) instead of 0(1/N) as it was the case for 
generic values of d. Let us write each symbol entering Eq. fliTD as 

#=*-+#+°G0. (36) 

and expand the /3-function in powers of 1/N to 0(1/N). Then comparing with Eqs. ( |T5D 
and (|29|), one finds the relations 

1 

a 2,oo + 7l,oo + Ol, 00 + Cl,oo — 

5 1>00 {A 2>1 - 2Ai) + Ci,oo(A 3 ,i - 2Ai) - 7i t00 Ai = B(2). (37) 
The leading exponent A is known to 0(1/N 2 ) for d = 3 [4~I,4~2 



A " 1 " T*N 3^ TP + ° [w) • (38) 

Eqs. (p7|) show that at least one among the coefficients 71, 5\ and Ci must be non-zero, 
and therefore that 0:2,00, 7i,oo, ^i,oo, and Ci,oo are discontinuous at d = 3. This could be a 
feature of the large- N limit: for finite values of N it is still possible that all coefficients be 
continuous in d f43|j . 

In conclusion our explicit calculation shows that non-analytic terms are present in the j3- 
function defined in Eq. (|T3J). Notice that this result is valid for all dimensions with 2 < d < 4 
and thus also in the e-expansion when one uses a massive renormalization scheme. In this 
case the singularity is of the form (g* — g) 2//e , a behaviour which has also been predicted 
from a large-order analysis of perturbation theory |44| , f45|| . However, for the e-expansion, 



it is still possible that, as conjectured in Ref. |j4q| , the /3-function is analytic in a massless 
renormalization scheme. The question requires further investigation. 



For small values of N, i.e. N = 0, 1, 2, 3, the renormalization group analysis of Ref. f47 
shows that A 2 /A ~ 2, A 3 /A ~ 3, and 1 + 1/A ~ 3. The closeness of such values to 
integer numbers may explain the small effects of the confluent non-analytic corrections in 
the procedures used to estimate g*. 
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IV. e-EXPANSION RESULTS 



A. Computation of g* to order 0(e 4 ) 

In this Section we will give the explicit expression of g* up to four loops, i.e. 0(e 4 ). 
There are essentially two different methods to perform the calculation. One may follow the 
previous section: starting from the definition, one computes the renormalized two-point and 
four-point functions, then derives the /3-function and finally obtains g* from the equation 
(3{g*) = 0. Alternatively, one may compute g at three loops in terms of g^g and e. Then, 
to express g* in terms of e only, one can use the four-loop expression (i.e. 0(e 4 )) for the 
fixed point value of g^r [fj8| . We have followed this strategy, which allows us to obtain g* to 
order 0(e 4 ) calculating only three-loop graphs. Some intermediate results are presented in 
App. ||. In the framework of the e-expansion, we found convenient to consider the rescaled 
coupling g, defined by Eq. (|6|). Due to the 0(e _1 ) factor multiplying g in the definition of 
g, the 0(e 4 ) of g* corresponds to the 0(e 3 ) of g*. The e-expansion of g* to 0(e 3 ) is given by 

g*(e) = J29ne n (39) 

?1=0 

with 

90 = 1, (40) 
_ 3(3A + 14) 

91 = (N + 8) 2 ' (41) 

-2864.85 - 1086.88A - 119.599A 2 - 7.07847A 3 
h = ^p^j , (42) 

298347 + 165854A + 38100.1A 2 + 4733.16A 3 + 240.959A 4 + 3.31144A 5 . . 

h = (nTW • (43) 

We have checked, see App. 0, that these expressions reproduce the 0(1/N) results of 
Ref. p|. 



B. Analysis of the e-expansion of g* 

The purpose of the calculation is, of course, the determination of g* in three dimensions 
and, possibly, in two dimensions. It is well-known that one does not obtain reliable estimates 
by simply setting e = 1 (or e = 2) in the corresponding series (0), since the expansion is 
strongly diverging. 

We have analyzed the e-series using the methods proposed in Ref. ||. The resummation 
technique is based on the knowledge of the large-order behaviour of the series. It is indeed 
known that the coefficients g n of the series g*(e) behave as 

g n = c(-a) n T(n + b + 1) (1 + 0(l/n)) . (44) 

The constant a, which characterizes the singularity of the Borel transform of g*(e) does not 
depend on the specific observable; it is given by [ P^ , pTH 
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3 

" N + 8' 



(45) 




(46) 



fc=0 



we have generated new series R p (a, b; e) according to 



R p (a,b;e) = ^B k (a,b) I dte~H h 



u{et) k 



(47) 



fc=0 



(l-u(et)) 



a 



where 



u(x) 



a/1 + ax — 1 



(48) 



\/l + o,x + 1 



The coefficients B k (a,b) are determined by the requirement that the expansion in e of 
R p (a,b;e) coincides with the series (|46"|). For each a, b and p an estimate of R is simply 
given by R p (a, b;e = 1). 

For g* we have computed the series ([17]) for many values of a and b, and for p = 2 and 
p = 3, obtaining in this way many different estimates of g* in three dimensions. We have 
noticed that, for a > 1, the estimates strongly oscillate with the number p of terms which 
are considered. These oscillations increase in size as a increases. For this reason we have 
decided to keep a in the interval — 1 < a < 1. Then for each value of a and iV we have 
considered various choices of b. In each case we have found an integer value of b, b opt , such 
that 



In other words b opt is the integer value of b that minimizes the difference between the 
estimates from the 0(e 2 ) and 0(e 3 ) series. In a somewhat arbitrary way we have then 
considered as our final estimate the average of R p (a,b;e = 1) with — 1 < a < 1 and 



Of course the real problem is the determination of the error bar. We have decided here 
to use an algorithmic procedure. The reason is that, if error and mean value are determined 
algorithmically, there is less chance to introduce unwillingly a systematic bias due to what 
we expect to be the "correct" value. Of course we do not only want to determine the 
error bars in a completely automatic fashion, we also want to have error bars which are 
reasonable. Since, as we will discuss below, we will obtain many different estimates of g* 
from the e-expansion, the basic requirement will be that all estimates should be compatible 
among each other. Notice that this is a requirement of internal consistency only, which 
does not use any external information. In this sense our procedure will be totally unbiased. 
Therefore discrepancies with results obtained from other methods will be meaningful. We 
have thus fixed our error bar as the sum of two terms: the first one is the variance of the 



#3(a, b op t] e = 1) w R 2 (a, b opt ; e = 1). 



(49) 



2 + b opt <b<2 + b, 



'opt- 
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values of ^(a, b; e = 1) with — 1 < a < 1 and [b opt /2 — lj < b < \3b opt /2 + 1] ; the second 
is the difference between the estimates from the series at order 0(e 3 ) and 0(e 2 ). 

To understand the reliability of our method we have reanalyzed the series of the critical 
exponents to order 0(e 5 ) |48| , [5T1J52]| and compared our estimates and error bars with the 
results of Refs. [ p3| , |54| . For N — 1, by analyzing the series of Ref. [pT]| , the estimates 
v = 0.6305(25) and 7 = 1.239(4) were obtained in Refs. |53|j54| . Analyzing the same series 
by our method we find v = 0.6272(32) and 7 = 1.2343(35). These estimates agree within 
errors and also the error bars are very similar. Our algorithmic procedure appears to give 
results similar to the estimates of other authors^. 

It has been noticed that the estimates of the e-expansion can be improved if one knows 
the (exact) value of the quantity one is considering in two dimensions [15[| or even in one 
dimension [p5] . We will now try to do more. First of all we expect g*(e) to be analytic in 
the domain < e < 4. This conjecture can be checked in the large- N limit using the exact 
results of Ref. |T2j . Moreover it has been implicitly assumed in the dimensional expansion 
around d = done in Refs. 0,0. Then we will try to improve our estimates using the 
exactly known results for e = 3 (d — 1) and e = A (d = 0) and furthermore the estimates for 
e = 2 (d = 2). 

Let us now discuss the procedure which is a simple generalization of the technique pre- 
sented in Ref. fll5| . Suppose that for e = e\ the exact value R ex (ei) is known. One may then 
define 

'R(e)-R^(e 1 )' 



R(e) 



(50) 



and a new quantity 



^imp(e) = iUei) + (e - ei)R{e). (51) 



New estimates of R at e = 1 can then be obtained by applying the resummation procedure 
we described above to R(e) and then computing i? imp (l). This strategy can be generalized 



4 It should be noticed that the series of Ref. [51] contained an error in the five-loop coefficients 
as shown in Ref. [52]. An analysis of the correct series has never been published. In Ref. [52] 
the authors only mention that the analysis of the correct series gives results that are consistent 
with those of Refs. []53,^4|. For the sake of completeness, we performed a new analysis of these 
series using our procedure. For N = 0, 1 we employed a constrained analysis using the known 
exact results in two dimensions, as in Refs. |53|,54|. We found v = 0.5882(11) and 7 = 1.1559(10) 
for N = 0; v = 0.631(3) and 7 = 1.240(5) for N = 1. We mention that perfectly consistent 
results are obtained using also the homografic transformation e' = Ae/(A — e) || with A = 4,3 
for N = 0, 1 respectively. We also report results from an unconstrained analysis of the series for 
N = 2,3: v = 0.664(3) and 7 = 1.304(7) for N = 2; u = 0.699(4) and 7 = 1.372(6) for N = 3. 
For iV = 2, 3 this analysis seems to underestimate (slightly) the values of the exponents (see e.g. 
Refs. P,|60||). Indeed by employing a homografic transformation with A = 2 (or constraining the 
exact two-dimensional values: \jv = I/7 = 0) one obtains larger values by approximately two per 
cent {y ~ 0.676 and 7 ~ 1.324 for N = 2, and v ~ 0.712 and 7 ~ 1.394 for iV = 3), but with much 
larger "errors". 
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to an arbitrary number of points: if exact values R ex (ei) 
dimensions E\, . . ., k > 2, then one defines 



Rex( e k) are known for a set of 



Q(e) = E 



2 = 1 



(r - o) . n 



(52) 



and 



and finally 



(53) 



-Rimp(e) 



g( e ) + J R( e )ln(e-e 



(54) 



i=l 



One can easily verify that the expression 



Q(e) + J R(0)ln 



e - €i 



(55) 



i=l 



represents the fc-order polynomial interpolation among the points e = 0, ei, e&. Again 
the resummation procedure is applied to R(e) and the final estimate is obtained computing 
-Rimp(l)- The idea behind this method is very simple. If, for instance, the value of R for 
e = 2 is known, one uses as a zeroth order approximation at e = 1 the value of the linear 
interpolation between e = and e = 2 and then uses the series in e to compute the deviations. 
If the interpolation is a good approximation one should find that the series which gives the 
deviations has smaller coefficients than the original one. Consequently also the errors in the 
resummation are reduced. In our case the value of g* is known for e = 3 and e = 4. For 
e = 2 we will use estimates obtained using strong-coupling methods or, for larger values of 
N, derived from the 1/N expansion. Using this additional information we will be able to 
substantially reduce the error on our estimates. 

As an example of the procedure let us consider the Ising model, N = 1. 

In this case the series of g*(e) is given by 



g*(e) = 1 + 0.629629e - 0.621613e 2 + 0.954535e 3 + 0(e 4 



(56) 



The coefficients are big, and at first sight it may seem hopeless to try to get an estimate of 
g* for e = 1. However the series alternates in sign and therefore one may obtain reasonable 
results after a Borel transformation. Indeed we find g* = 1.37 ± 0.09 which is already quite 
good. We can now try to improve our estimates using the value of g* for e = 2. Using 
the two-dimensional estimate g* = 1.7540(2) (see next Section) and Eq. (|54|), we get a new 
series for g*: 



'imp 



1.7540 + Ag(l), 



(57) 



where 
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Ag(e) = 



-0.37698 + 0.12632e - 0.24765e 2 + 0.35345e 3 + 0(e 4 ). 



(58) 



The coefficients of the new series are on average a factor of three smaller than the coefficients 
of the original series and the simple interpolation already gives a good estimate of g* mp : 
indeed 1.754 + Ag(0) = 1.377, not very far from the correct value. Consequently one 
expects a corresponding gain in the error bar. Indeed we get 1.400 ± 0.017. There is also 
an additional error due to the uncertainty in the two-dimensional result which in this case 
however turns out to be negligible. One can go further and use both the estimate for e = 2 
and the exact result for e = 3 [12|, which is g* = 9/4. Following the procedure we presented 



above, we get a new estimate of g* from 



'imp 



1.2580 + Ag(l). 



where 



Ag(e) = 0.07937 + 0. 110671c - 0.128206e 2 + 0.192895e 3 + 0(e 4 



(59) 



(60) 



On average the coefficients of the new series are a factor of two smaller than those of the 
series in which only the estimate at e = 2 was used. Correspondingly we obtain a slightly 
more precise estimate g* = 1.395 ± 0.016. Again the error due to the uncertainty on the 
two-dimensional result is negligible. 

Finally we can include also the known value of g* for e = 4, g* = 3. Using the values at 
e = 2, 3, 4, we get a new expansion in the form 
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imp 



1.514 + Ag(l). 



where 



Ag(e) = -0.13095 + 0.05027e - 0.08359e 2 + 0.12377e 3 + 0(e 4 



(61) 



(62) 



Apart form the first term of the series, all the other terms have coefficients which are smaller 
than those of the series (^) and thus we expect an additional reduction of the error. Indeed 
we find g* = 1.397 ±0.008. 

One can try other possibilities, constraining the series only at e = 3 or e = 4, or at 
any possible pair. In Table [I] we present also the results obtained by constraining the series 
at e = 3 and at e = 3,4. The estimates are all consistent among each other, thus giving 
confidence to the final result. 

We have applied this method to all values of N. In one dimension g* was computed in 
Ref. [0 finding 
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and 



N + 2\ AN) 

_* N + 8 



for JV>1, 



for N < 1. 



(63) 



(64) 



In zero dimensions we generalize the result of Ref. |57| to any integer N > 1. We get 
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(65) 



It is not clear how to obtain the value of g* for N = 0. The one-dimensional results, Eqs. fl63|) 
and (0), indicate that one cannot naively set N = in the formula obtained for N > 1. The 
incorrectness of this analytic continuation can also be verified by our constrained analysis: if 
one uses the prediction g* = 4 for e = 4 one obtains results which are in total disagreement 
with the other estimates. 

In d = 2 we will use the best available estimates. For iV < 4 we consider the estimates 
obtained from the analysis of the strong-coupling series, more precisely, g* = 1.679(3), 
1.7540(2), 1.810(10), 1.724(9), 1.655(16) respectively for N = 0,1,2,3,4 (see next Section). 
For iV > 8 we use the large- N expression [|l^] 

JV + 8 / _ O602033N 
y iV + 2 V N J K ; 

Of course the prefactor in front of the previous expression is arbitrary. However if one uses 
this particular form one finds a small 1/N correction and good agreement up to iV = 4. 
Indeed comparison with the (imprecise) strong-coupling result at iV = 4 we have reported 
above shows that Eq. (|66D reproduces the correct result with an error smaller than 4%. 



In the absence of better estimates for iV > 4 we have thus used formula (36) with an 
"estimated" error of 0.64(?*/iV 2 (which reproduces the difference of 4% found at N = 4). 

The final results of our analysis for selected values of iV are reported in Table [IV]. 

Our results can be checked in the large- iV limit using the exact result, ~g* = 1 + 4.4540/iV 



T2||. Let us consider the constrained analysis in d — 0, 1, 2, that is the one that provides the 



most precise estimates. In this case g* is determined from 



where 



, 4.9439 . 

1 + — ^ + A<?(1), (67) 



Ag(e) = —(-1.09695 + 0.90986e - 0.359826e 2 + 0.051144e 3 ). (68) 



Applying our resummation procedure to Ag(e) we obtain the estimate 

S-l + iffl (69) 

in perfect agreement with the exact result reported above. 

We have also repeated the analysis in two dimensions. In this case of course it is more 
difficult to get precise estimates: the unconstrained expansion gives results with errors 
of order 20 — 50% and it is therefore practically useless. Better estimates are obtained 
constraining the expansion in one and zero dimensions. The results for these two cases 
are reported in Table O, and are consistent. Of course the errors are larger than in three 
dimensions, but still small if one considers the series we started from (for = 1 see Eq. (|56|) ). 
The resulting estimates of g* are in good agreement with the strong-coupling and large- N 
estimates we used above in the analysis of the three-dimensional case, thus supporting their 
use. 



15 



V. STRONG-COUPLING EXPANSION 



In order to understand the reliability of our previous estimates for the renormalized 
coupling constant, it is useful to compare them with the results of a quite different approach, 
such as the analysis of the high-temperature expansion for the iV-vector model. In this 
context it is convenient to consider the coupling / defined in Eqs. (§) and (]XT|) . The finiteness 
of /* is related to the hyperscaling relations of the critical exponents. A detailed analysis 
of the strong-coupling expansion for all values of N was presented in Ref. [E3] (see also 
Ref. pip extending the numerous studies for the Ising model |T7|- pO| , p2| , p3| . The results 



were in substantial agreement with the renormalization-group estimates, except for small 
values of N, where relatively small systematic deviations were found. The reason of these 
discrepancies is the presence of confluent singularities at (3 C . Indeed, in general /(/?) behaves 
as 



f(P) = f* + c A ((3 c -(3) A + 



(70) 



close to the critical point. The traditional methods of analysis, like Pade and Dlog-Pade 
approximants, are unable to handle an asymptotic behaviour like ([70]) when A is not an 
integer number, thus leading to a systematic error. 

To take into account this kind of confluent corrections in full generality, one should 
consider integral approximants, which, however, require long series to detect non-leading 
effects, and in practice need to be biased to work well. Roskies proposed a rather simple 
method to handle the leading confluent singularity in the Ising model where A ~ 1/2. 
He showed that the effect of the non-analytic terms can be significantly reduced by a suitable 
change of variables. Equivalently one can use suitably biased integral approximants | 60fl . In 
Ref. [2j| the series for the renormalized coupling constant for the Ising model was analyzed 
biasing A = 1/2. The estimated value of /* was significantly lower than previous estimates, 
and now in good agreement with the renormalization-group prediction. This procedure was 
also successfully applied to the calculation of the critical exponent of the specific heat from 
the low-temperature expansion [p8|], providing results consistent with field theory, while a 



standard analysis neglecting confluent singularities led to a quite inconsistent estimate [59 



We have decided to repeat the analysis of Ref. |I2j for all values of N using similar ideas. 
The new results show systematic differences from the old ones and now they are in much 
better agreement with the renormalization-group estimates. 

The strong-coupling expansion of f{(3) has the following form 



In three dimensions the available strong-coupling series allow us to calculate 

A{fi) = P 3/2 f(P) 



(71) 



(72) 



up to 14 th order (using x and m 2 up to 15th order 0,^T|], and %4 up to 14th order |)2]j63|D . 
Longer series are available for the Ising model, N = 1. On the cubic lattice, using the 
published series [Q.[l5|.|3|] , one can derive A((3) up to 16th order. Moreover series on 
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other lattices are available, which allow us to calculate A(J3) on the b.c.c, f.c.c. and di- 
amond lattice, up to 13th, 10th, and 19th order respectively (using series published in 
Refs. pi 1 53 , 0] , 56[ , pl| ) . We mention that longer series for all values of N on the cubic and 



b.c.c. lattice have been announced in Refs. |2T[|60ll , but they have not been published yet. 



The idea of the Roskies transform (RT) |32[ is to perform biased analyses which take 
into account the leading confluent singularity. For the Ising model, whereQ A ~ 1/2, one 
replaces the variable (3 in the original expansion with a new variable z, defined by 

1 - z = (1 - /V/? c ) 1/2 . (73) 

Of course a quite precise estimate of (3 C is required here. If the original series has square-root 
correction terms, the transformed series has analytic correction terms, which can be handled 
by standard Pade or Dlog-Pade approximants. 

In order to analyze models with different values of N and therefore with A ^ 1/2, one 
can generalize the Roskies transformation and consider the change of variable 



1- Z = (1-P/(3 C )». (74) 

In the following we will refer to this trasformation as GRT. For b = A, this mapping makes 
the first correction to scaling analytic. Thus standard approximants can handle it correctly. 
Of course non-analytic terms still survive due to subleading corrections, but they should be 
less important as far as A is sufficiently smaller than the next exponents. Anyway they still 
represent a (hopefully small) source of systematic error for our analysis. From Eq. ( [74] ) it 
follows that 

f{fi) — f(z) = z-^A{z). (75) 

In order to estimate f(/3 c ) = f(l) we analyze the expansion of A(z) in powers of z. As we 
shall see, the use of the mapping ([74]) leads to a much better agreement with field-theoretic 
estimates. 

Note that the relevant singularity of f(z), i.e. the one at z — 1, is no longer the closest to 
the origin. Indeed, the antiferromagnetic singularity at — f3 c is mapped closer to the origin 
than that at j3 c , but still at negative values of z. We expect its effect to be small when 
evaluating the resummed series around z = 1. 

By employing standard resummation methods, we studied the behavior of f(z) around 
z — 1. For the sake of comparison, we also performed a standard analysis, i.e. without 
using the mapping ([74]). We constructed various types of approximants to the series of 
A(z), such as Pade approximants (PA's), Dlog-Pade approximants (DPA's) and first-order 
inhomogeneous integral approximants (IA's) (for a review on the resummation techniques 



see for example Ref. |70|). Note that in principle IA's should be able to detect the first non- 



analytic correction to scaling in Eq. ([70]) , but they probably need more terms of the series, 



5 We mention the recent estimate of A for N = 1 obtained by the fixed-dimension field-theoretic 
method: A = 0.498(8) |37|. Notice that this is somewhat lower than the high-temperature esti- 
mate MM A = 0.54(3). 
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and practically need to be explicitly biased as in the case of PA's and DPA's. Indeed the 
IA results without the GRT turn out to be substantially equivalent to those obtained from 
PA's and DPA's. In all cases we considered only quasi-diagonal approximants^. Then we 
evaluated the approximants at z = 1 in order to obtain an estimate of the fixed-point value 
of /. Very precise values of /3 C , which are needed for the mapping (|73|), are available in the 
literature from different calculations (see for example Refs. [|7T H75| , |60| ,p^1). Errors due to the 
uncertainty on the value of (3 C turned out to be negligible in our analysis. As estimates of A 
we used the field-theoretic prediction for iV < 24, and its large- N expression fl38|) for larger 
values of N. Since at N = there is a relatively large difference between the field-theoretic 
prediction, A = 0.470(25), and the Monte Carlo estimate, A = 0.515(7) |7§, in our 
analysis we considered the very conservative value A = 0.50(5). Similarly for N = 1 (see 
footnote |5|) we consider the value A = 0.50(5). 

In a PA or DPA analysis of quantities with a confluent non-analytic correction, the 
singularity should be mimicked by shifted poles at (3 > (3 C . PA's and DPA's of A((3) present 
indeed singularities typically at (3 ~ 1.1 4- 1.2 /3 C . On the other hand, the approximants of 
the series in z (cf. Eq. ([74])) do not show singularities close to the new critical value z — 1, 
conferming the effectiveness of this change of variable. Most approximants of A(z) do not 
present singularities near the real axis in the region Re z < 1.5. 

Table [III] shows the results of our strong-coupling analysis]] for selected values of N. 
There we report also the values of f3 c and A we used. The results of the analyses which 
use the GRT are quoted with two errors: the first one is the spread of the approximants for 
b = A, the second one is due to the uncertainty on A, and it is obtained by varying b. We 
note that for N = 0, 1, 2 the estimate of /* increases with increasing b in the GRT. 

A few comments on the results of Tables III are in order. 

(a) For most values of N, the difference between the results of the analysis with and 
without the use of the GRT, although relatively small, is larger than the apparent error of the 
single analysis (which is estimated by looking at the stability of the different approximants). 



6 Given a nth order series, we considered the following quasi-diagonal approximants: [l/m] PA's 
and DPA's with I + m > n — 2 and l,m> § — 2 (I, m are the orders of the polynomials respectively 
in the numerator and denominator of the PA of the series in the case of PA's, or of its logarithmic 
derivative in the case of DPA's); [m/l/k] IA's with m + I + k + 2 = n and [(n — 2)/3j — 1 < 
m,l,k < \(n — 2)/3] +1 (m,l,k are the orders of the polynomial Q m , Pi and Rk defined by the 
first-order linear differential equation Q m {x)f'{x) + Pi{x)f{x) + Rk{x) = O ^x k+l+m+2 ^j , whose 
solution provides an approximant of the series at hand) . 

7 As estimate of /* from each class of approximants (i.e. PA's, DPA's, and IA's) we took the 
average of the values at z = 1 (or (3 = (3 C when performing a standard analysis) of the non-defective 
approximants using all the available terms of the series. The error we quote is the square root of 
the variance around the estimate of the results from all the non-defective approximants listed in 
the footnote |6|. Approximants in the generic variable x are considered defective when they present 
spurious singularities close to the real axis for Re x < x c . The results from PA's, DPA's, and IA's 



are then combined leading to the estimates shown in Table III 
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This indicates that the systematic error due to the neglecting of confluent singularities is 
much larger than the error obtained by a stability analysis of the results. 

(b) The estimates from the GRT analyses are globally in much better agreement with 
the field-theoretic predictions than the results obtained neglecting the confluent singularities 
(except for N = 4, for which the agreement of the non-biased result was already satisfactory). 
This can be seen in Table [IV] where the corresponding results for g* are reported. We 
note that the strong-coupling estimates are systematically slightly higher for N > 3. The 
uncertainty of the GRT results is approximately 1%, thus providing an accurate check of 
the field-theoretic calculations by a different approach. 

(c) For the Ising model, universality among formulations on the cubic, b.c.c, f.c.c, and 
diamond lattices is well verified by the results of our GRT analysis. Assuming universality, 
our overall estimate is /* = 23.55(15). This is consistent with the result of the biased 
analysis (sligthly different from ours) of the strong-coupling expansion on the cubic lattice 
of Ref. [^j, /* = 23.69(10). We note that universality is apparently shown also by the 
results of the standard analysis, although they lead to a different estimate of /*. This may 
be explained by noting that if the values of the non-universal coefficients ca in Eq. ( |7Tj ) are 
approximately the same for all the considered lattice formulations, they may give rise to 
similar systematic errors leading to an apparent universality of the results. This would not 
be surprising. Indeed, already the leading amplitudes of many non-universal quantities have 
close values in various nearest-neighbor lattice formulations (see e.g. the results reported 



in Ref. [22], and the estimates of some amplitudes of leading scaling corrections reported 



in Ref. ||77|| ). In Ref. f22| a different analysis still neglecting confluent singularities (where 



the amplitudes of Xi m 2 and %4 were indipendently calculated by a first order integral 
approximant analysis to give an estimate to /*) led to the following results: /* = 24.55to;<25 
on the cubic lattice, /* = 24.39(9) on the b.c.c. lattice, and /* = 24.50(13) on the f.c.c. 
lattice. Universality is nicely observed, but, again, the final result is larger than the GRT 
estimate, and therefore also than the field-theoretic ones. Notice that these estimates of /* 
are lower than our results obtained without using the GRT. This is probably due to the 
different procedure used to estimate /*, and to the fact that different series were analyzed, 
(d) With respect to the standard analysis, the GRT results are in much better agreement 



with the formula obtained by a 1/N expansion [12 

/* = 16tt 



1.54601 / 1 



(76) 



This equation gives: /* = 48.646 for N = 48, /* = 47.837 for N = 32, /* = 47.027 for 
N = 24, /* = 45.408 for N = 16, etc... This agreement provides additional support to 



the formal argument presented in Sec. |T|, according to which analytic terms should not be 
present in the expansion of / around (3 C . Indeed, in the opposite case, one would not expect 
a substantial improvement using the mapping ([74] ) for iV sufficiently large, since A — > 1 for 
N -> oo. 

We mention that for the Ising model high-temperature techniques have also been used 
to obtain a dimensional expansion of the Green's functions around d = 0. The analysis of 
these series presented in Ref. |5?| led to the quite good estimate /* = 23.66(24). 

Let us compare our strong-coupling predictions with the results obtained from Monte 
Carlo simulations. Monte Carlo estimates of /* for the Ising model on the cubic lattice can be 
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found in Refs. f24}-|31~|| . In Fig. |] we compare some of the data of the most recent works with 
the Pade resummations of the strong-coupling series with and without the GRT. The data of 
Ref . pij , which are those closest to criticality, have been obtained by employing a finite-size- 



scaling technique, which allowed the authors to get data up to a value of (3 corresponding 
to £ ~ 30. These data show f(j3) apparently flattened around 24.5(2). Ref. |3(| presents 
data up to £ ~ 10, from which the authors obtain the value /* = 25.0(5). Fig. [I] shows 
that Monte Carlo data are in substantial agreement with our GRT analysis^. But biased 
approximants extrapolate to a smaller value of /*. It is worth mentioning that the use of 
the GRT to bias the strong-coupling approximants is, in a sense, equivalent to the use of 
the function 

/(/?) = r + c A (/3 c -/?) 1/2 (77) 



for the extrapolation to j3 c of the Monte Carlo data at /3 < f3 c . Using the function (|77|) to 
fit the Monte Carlo data of Fig. Q one gets /* = 23.7(2) and c A = 34(3) with x 2 /d.o.f ~ 0.6 
(here we assumed all data to be independent)^ which is perfectly consistent with the strong- 
coupling and field-theoretic estimate of /*. Finally we mention the result of Ref. PU 



/* = 23.3(5), obtained by studying the probability distribution of the average magnetization. 

We mention another application of the Roskies transform Q73D , that is the analysis of the 
low-temperature expansion of the quantity u defined in the broken phase of the Ising model 

by 

(78) 

where M is the magnetization, u plays the important role of a zero-momentum low- 
temperature renormalized coupling constant in the study of the 4 theory directly in 
d = 3 [f78|| . The most precise determinations up to now have been apparently obtained by 
Monte Carlo simulations ]79| (where data have been fitted by using Eq. (|77"D ): u* = 14.3(1), 
and by two different analyses of the low-temperature expansion ||50| , [22"1 , which lead to appar- 
ently inconsistent results: u* = 14.73(14) f8(J and u* = 14.14(14) p2jpl| . We repeated the 



analysis of the low-temperature expansion using the RT. The series published in Refs. 
allow us to calculate the expansion of u in powers of e _4/3 up to 21th order. Quasi- diagonal 
PA's of the Roskies-transformed series give u* = 14.3(1) (without using the RT one obtains 



u* = 14.7(1) as in Ref. ||80||), thus confirming the Monte Carlo result of Ref. [79| and the 
low-temperature analysis of Ref. |^2| . 

In Sec. [IV] we used estimates of /* in two dimensions for our constrained analysis of the 
e-expansion. In the following we shortly discuss their derivations. An analysis with Pade 



approximants of a 17th-order series was presented by Butera and Comi MM. For N <2 they 



8 One should also take into account that the data of Ref. [31| at different /?'s are not statistically 
independent because they have been obtained by a finite-size scaling technique. 

9 If we do not include the data for the lowest value of /?, corresponding to £ ~ 3, we obtain 
/* = 23.9(3) and c A = 27(7). 
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found: /* = 10.53(2) for N = 0, f* = 14.693(4) for N = 1, and /* = 18.3(2) for N = 2. 
We reanalyzed the same series using also DPA's and IA's. For N = 2 we also considered 
the series in the internal energy E [O]. We found results in total agreement: 10.55(2) for 
N = 0, 14.693(1) for N = 1 and 18.2(1) for N = 2. Moreover, using the series for the Ising 
model on the triangular lattice published in Refs. p4] , |65| , |85| , we obtained the strong coupling 
expansion of (3f((3) to 14th order. Its analysis gave /* = 14.695(1). A comparison with the 
square-lattice result (i.e. assuming universality) leads to the final estimate /* = 14.694(2) 
for the two-dimensional Ising model. 

For N = 3 no estimate was reported in Ref. pi|. Our analysis gives 



/* = 19.7(1) for N = 3. (79) 

For comparison we mention the field-theoretic estimate obtained in Ref. [|f| for N = 3: 
f* = 20.0(2). 



Estimates for several values of iV > 3 were reported in Ref. [j34 |. For N = A f* = 20.9(1) 
was found, which is in good agreement with our reanalysis: /* = 20.8(2). For larger values 
of N — we will be interested in N > 8 — the analysis of the strong-coupling series gives 



results with a somewhat large error and in this case the 1/N expression |T2] 

f* = 8n 



1 0.602033 ^ / 1 



(80) 



should provide more precise estimates. Indeed it is already a good approximation for N = 3 
and N = 4, where it gives /* = 20.09 and /* = 21.35 respectively. 

The reader should notice the small uncertainty of the estimates for N = 3, 4 in spite of 
the fact that (3 C = +oo. This is due to the fact that, according to field theory, f(/3) like any 
dimensionless renormalization-group invariant quantity behaves as 

m-r~^, pi) 

for sufficiently large f3. Hence the corrections to /* decrease exponentially in (3. This 
important point was overlooked in Ref. ||84|| . As a consequence of Eq. ( [Sip the scaling 



region, where the function f(/3) approximately reaches the asymptotic value, may begin 
quite early. One may obtain good estimates of the dimensionless renormalization-group 
invariant quantities already at £ > 10, which is still within the reach of the strong-coupling 
extrapolation, at least for not too large values of N, say N = 3,4 [|S^| . For instance in 
Fig. |2| we show f(/3) for iV = 3, 4 versus the correlation length. Scaling is nicely verified by 
our strong-coupling calculations. Our estimates of /* at N = 3,4 are obtained from f(/3) 
at £ « 10. The behaviour ( |31~D explains the success of the strong-coupling method when 
applied to dimensionless renormalization-group invariant quantities (for other examples see 
Ref. |U). 



We finally mention the estimates of /* for N = 2 and N = 3 obtained by a Monte 
Carlo simulation together with a finite-size scaling extrapolation f89 |. For N = 2, fitting 



the data of Ref. with £ > 10 to a constant, we get /* = 17.7(3). For N = 3 the result is 
/* = 19.8(3). They are in good agreement with the strong-coupling results presented above. 
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VI. CONCLUSIONS 



We have studied some issues concerning the fixed-point value of the zero-momentum 
four-point renormalized coupling g in O(N) models. The coupling g plays an important role 
in the field-theoretic perturbative expansion at fixed dimension, which provides an accurate 
description of the symmetric phase. In this approach the value of g* is essential to compute 
the critical exponents, which are obtained by evaluating appropriate anomalous dimensions 
(calculated as functions of g) at g*. 

The first important issue we have discussed is related to the presence of confluent singu- 
larities at the zero of the Callan-Symanzik /^-function. In order to understand this problem 
we have considered the framework of the 1/N expansion, which provides an analytic and non- 
perturbative approach. The analysis of the next-to-leading order of the /9-function shows the 
presence of confluent singularities at its zero, as argued by Nickel ||11|| . In generic dimensions 
d 7^ 4 — 2/n (with n G IN), the leading non-analytic corrections are 0(1/N) and are related 
to the exponents A 2i 3/A and/or 1 + 1/A. Since they are degenerate at N = oo, one cannot 
distinguish them by an 0(1/ N) calculation. No term associated with an exponent 1/A is 
found. In three dimensions one meets a phenomenon analogous to the so-called Abe-Hikami 
anomaly []13 |. One indeed finds that for d = 3 the non-analytic contributions are 0(1) in the 



1/N expansion, even if the /3-function appears analytic to leading order. This is essentially 
due to a further degeneracy occurring for d = 3 at N = oo among A 2 3 /A, 1 + 1/A and the 
analytic correction with exponent two. 

In the analysis of the ^-expansion performed by Le Guillou and Zinn- Justin || with 
an additional hypothesis of analyticity, the presence of such singularities may cause a slow 
convergence to the correct fixed-point value, thus leading to an underestimation of the real 



uncertainty An accurate check of the (^-expansion results was our major motivation for 
the extension of the e-expansion of g* to 0(e 4 ), and for a reanalysis of the strong-coupling 
expansion in the lattice iV-vector models. 

We obtained rather accurate estimates (with an apparent precision of approximately one 
per cent, see Table [TVD from the analysis of the 4th order e-expansion of g* that exploits the 
known values for 0(N) models in lower dimensions, g* is indeed expected to be analytic in 
the domain < d < 4, as can be verified in the large- N expansion to 0(1/ N). We plan to 
extend this analysis to the low-magnetization expansion of the effective potential, which is 
parametrized by the zero-momentum n-point renormalized couplings. 

The agreement with the ^-expansion estimates is globally good. For N > 2 there is 
full agreement. The results for N = 0,1 are slightly lower than the estimates given by Le 
Guillou and Zinn- Justin |J, thus favouring the more general analysis done by Nickel |§. 
This would lead to a small change in the estimates of the critical exponents, since they 
depend crucially on the value of g*. For instance, consider the case N = 0. For this model 
a very precise estimate^ of the exponent 7 has been recently obtained by a Monte Carlo 
simulation: 7 = 1.1575(6) [pO| . On the other hand, the analysis of the ^-expansion, i.e. the 



10 It is worth mentioning that recently a very precise estimate of v was obtained by Belohorec and 



Nickel [76 1 by a Monte Carlo simulation of the Domb- Joyce model: v = 0.58758(7). 
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/^-function to 0{g 7 ) and the function 7(3) to 0(g 6 ), of Ref. g] led to 7 = 1.1615(20). A 
more precise estimate is reported in Ref. 0,7 = 1.1607(12): it is obtained using the same 
resummation method but one additional order in the series of j(g) [§]. Reanalyzing the 
0(g 7 ) series of 7(g), using the same method of Ref. ||, the authors of Ref. 0] reported the 
estimate 

7 = 1.1616 + 0.11(0* - 1.421) ± 0.0004, (82) 

where g* is kept arbitrary. Our result for g*, i.e. g* ~ 1.39, thus suggests a lower value for 
7, 7 ~ 1.158, in substantial agreement with the results of the Monte Carlo simulations and 
with the analysis of the e-expansion (see footnote ^). Of course more precision and therefore 
longer series are necessary to be conclusive. We also mention that a recent analysis of the 
21st order strong-coupling expansion biasing A to the known approximate value has given 
7 = 1.1594(8) on the cubic lattice, and 7 = 1.1582(8) on the b.c.c. lattice (slightly 
larger values have been obtained by unbiased analyses). 

As a by product of the e-expansion we obtained rather accurate estimates for the two- 
dimensional models, which are in good agreement with other estimates from lattice ap- 
proaches (both strong-coupling expansion and Monte Carlo simulations) and 1/A^-expansion. 

Finally we reanalyzed the strong-coupling expansion of three-dimensional lattice N- 
vector models. In order to get accurate estimates of g*, we employed an analysis able 
to handle the leading confluent singularity. For this purpose we used a generalization |33| of 



the Roskies method JHJ consisting in a appropriate change of variable. Final results have 
an apparent precision of approximately one per cent. We found good agreement with the 
field-theoretic estimates. At A^ = 0, 1 the results seem to favor the lower values of g* of the 
e-expansion. 
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APPENDIX A: ASYMPTOTIC EXPANSION OF LARGE-iV INTEGRALS 

In this appendix we give a few technical details on the computations of Sec. [TTJ. 
The basic problem is the determination of the asymptotic expansion for g — > g*^ of 
integrals of the form 

'^-Iwf^Mw (A1) 

where f{u 2 ) is a rational function of u 2 and 6 = (g^ — g)/g, for 2 < d < 4. Using the 
definition of n(-u), cf. Eq. (|23|) , one can rewrite 

U/,e) = (!)"*„ jf ( A2) 
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where is defined in Eq. (p8|) and 



6(u) 



u 



1 + iJ 



2\ d/2-2 



u 



2' 2' 2' 4 + u 2 



(A3) 



We are interested in the large-w behaviour of S(u). Standard identities for the hypergeometric 
function give 



5(u) 



2\ A ' d 
u) 



c[l + A(u)} 



where c is defined in Eq. ([27]) and A(u) has the following asymptotic expansion 

A-d oo , 



fc=l 



fc=i 



(A4) 



(A5) 



the coefficients a& and can be easily computed for any value of d. From Eq. ( |A4[ ) one 
immediately sees that, for d > 2, 5(u) ~ u d_4 for large values of it. 

For d < 4 the singularities of / n (/) are due to the large-it domain. Indeed 5(u) goes to 
zero for u — > oo and thus one cannot perform a naive expansion in powers of O, since for k 
large enough l k (f,0) diverges. To compute the asymptotic expansion, let us consider 



R{t,u;Q) 



exp 



2 ^ 
tu 



[0 + t^ditu)} 

and its expansion for t — > oo which has the generic form 

^ tu y2k- P {i-d) 



(A6) 



R exp (t,u;Q) 
Then consider for 2 < d < A 



^ r h,k,p 



h,k,p>0 



[9 + c(2/it) 4 - d ] n 



+h' 



(A7) 



/(« 2 



(e + *(«)) 



- - -R exp (l, 6) exp 



-c [ - 

u 



A-d 



(A8) 



Because of the exponential factor, the second term is integrable for u — > 0. One can expand 
the integral ( |A8l ) in powers of G. Indeed the subtraction guarantees that each term is 
expressed in terms of a convergent integral^. Thus the integral (|A8| ) gives rise only to 
analytic contributions. The non-analytic terms can thus be computed from the subtracted 
term. Defining 



n The reader could rightly be worried by the presence of an infinite series in (|A8|). Indeed all 
formulae should be intended in a formal sense. More precisely the procedure is the following: given 
K, to obtain the expansion of I n {f) to order Q K , one should consider in R exp (t,u;&) only those 
terms which, for = 0, decrease for u — ► oo less than or as u - d - K {^ d ) , This truncated expansion 
should be used in ( |A8| ). 
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^c(-) , (A9) 
we see that we must compute integrals of the form 

J n, a = J dx (eT^F (A10) 

with a > 0. The computation is now trivial as 

J la = Q a e e T(a + l)r(-a,6), (All) 
(-l) n_1 d 11 - 1 

Jn > a ~ (n-l)!d0n-i Jl - a ' (A12) 
where T(—a, 0) is the incomplete T-function |JT). For —>■ we have 

r ( - a ,e) = r(-a)-e-«i;A^. (Ai3) 

APPENDIX B: e-EXPANSION CALCULATION 

In this Section we report the results of our calculation of g*(e) to order 0(e 4 ), e being 
defined as e = 4 — d. 

The computation requires the determinations of the massive two-point and four-point 
functions for p 2 — > to three loops. The most difficult graphs are those reported in Fig. |3|. 
The calculation is straightforward albeit long. We thus simply report the results: 
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9A 5tt 2 


12e 2 


48 




32 96 
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1 




(I TV 2 ) 


3? ~ 


3? + 




I 12 24 , 



i ^2 



1 vr 2 , c , 5 2 , 5 4 C(3) 



(Bl) 



24/ 8 16 4 4 6 

(6) = (a), (B2) 

-■ / ■ i 1 11A 7r 2 S\ So Sa 3S-7 5 t/ - s /t^„n 

( C )-^ + -(^-T + 7,)+^ + ^ + ^ + T + y- T + --^C(3), (B3) 



(d) = -^-^-^- ^ + ^ + ^-^-^5-^6-^7 + ^(3), (B4) 

-^ + ^3-^(3). (B5) 



Each result should be additionally multiplied by m 3e iVj where is defined in Eq. ([28] 
and m is the mass. The constants A and Si are defined by the following integrals: 



/°° , logt 1 , /1\ 27T 2 

A ^-( W^tTT = 3 (3) ~ ~ ~ 1 ' 17195 ' («!) 



/•oo 

51 = / p 3 dp \ogpd{p) 2 {L{p) - 2 log p) « 1.0207, (B7) 

5 2 = / p 3 dpd{p) 2 / cfe log 2 (p 2 x(l - x) + 1) - log 2 (p 2 x(l - x)) « -0.8619, (B8) 
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dx 



xKoixfKAxf - 



-x/2 



X 



X \ * 
IE + log - J 



0.45077, 



S 4 = / p 3 dp dip) 2 (Up) 2 - 4 log 2 p) « 5.9622, 
Jo 

/•CO 

S 5 = / p 3 dp\ogp d(pf(L(p) - 2) ^ 0.18604, 

5 6 = / p 3 dpd(pf / rfa; log 2 (p 2 x(l - ac) + 1) « 0.21105, 

jo Jo 

/•OO 

5 7 = / p 3 dpd{p) 3 {L{p) - 2) 2 » 0.19533, 

Jo 



(B9) 

(B10) 
(Bll) 
(B12) 
(B13) 



where ^e is the Euler constant, ^e ~ 0.577216, i^o(^) and ifi(x) are modified Bessel 
functions, and 



L(p) = £log 



i 



p 



2 ' 



(B14) 
(B15) 
(B16) 



The "Mercedes" graph (/) requires more sophisticated techniques. To compute it, we used 
the method of Kotikov [|92]| . Let us define 



I(M,m) 



d d p d d q d d r 
{2it) d {2n) d (2tt; 



d D(p, M)D(q, M)D(r, M)D(p - q, m)D(p - r, m)D(q - r, m) 

(B17) 



where D(p,M) is the massive free propagator with mass M. For m = this integral was 
exactly computed in Ref. |)IJ in all dimensions d. The expansion for e — > is given by 



7T 



/(M, 0) = M-^iVj I -C(3) - - + 0(6 

Using the strategy of Ref. |92| we end up with 

^ 2 f(m,M)-f(0,M) 



(B18) 



i(M, M) = 3M 



-3e 



dm 



m 2 (3M 2 — m 2 

where /(m, M) is the following quantity (finite for e — > 0) 

"if(m, M; q 2 )J(m,M; q 2 



+ /(lM))(L + |log^ + 0(e 2 ; 



(B19) 



/(m,M) 



-M 2 iV, 



/•OO 

/ g^cig 
Jo 



q 2 + M 2 

K(M, M; g 2 ) J(m, m; q 2 ) 1 if(M, M, g 2 ) 



log 



m 



g 2 + m. 2 327T 2 g 2 + m 2 M 2 

The functions J and if can be easily computed for e = from the integrals 



(B20) 



26 



K(m,M;q 2 ) = J S^D(p, M) 2 D(p + q,m), 



(B21) 



J(m, M; q 2 ) = -{mM)^l 2 {^)- d l 2 T (2 - f ) + J tJ^^ M)D(p + q,m). (B22) 
A numerical computation gives 

' ' ' (B23) 



I(M,M)=N 3 d M-^ i^-C(3) + H 
with H = -0.9825. 

Collecting everything together, we can compute the first four coefficients in the expansion 
(H) of g*. Explicitly 



£0 = 1, 



9i 
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3(3iV + 14) 
(iV + 8) 2 ' 
1 



1 



-2N 3 + 58N 2 + 520N + 1224 - -(13iV + 62) (JV + 8) 2 A 

3 
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(iV + 8) 4 

-12(5iV + 22)(iV + 8)C(3) 

8 (ArTsW ( 4iy5 ~ 99N * + 5404Ar3 + 57572iV 2 + 225312iV + 341312 



(B24) 
(B25) 

(B26) 



+ 



(iV + 8) 5 
1 1 



24 (iV + 8) 4 
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(23iV 3 - 209iV 2 - 2954iV - 6580) C(3) 



(lliV 3 - 4338iV 2 - 45312iV - 109472)A + 960(2iV 2 + 55iV + 186)C(5) 



+ 



(N + 8Y 



4 2 

^-(5N + 22) - ^-(19iV + 62) + 16(N - 1)S 3 
15 3 



+2(19iV + Q2)S 4 - -{17N 2 + A10N + 1328)5 7 - 8(5iV + 22)H 



1 



(N 



13 

24^ + 6S 2 + 13(N + 2)S 5 + — {N + 2)S 6 



(B27) 



Numerical values are reported in Sec. |IV A . 



We can use the results for N — > 00 of Ref. [JT2J to check our expression. For large values 
of N, g* is given by 



g* = i + l/(3- c /)2 d - 1 + 8- / f^— -d{u) 



4 ^ ) + 9 (H^ 



(B28) 



where g*^ and U.(u) are defined in Eqs. ([H]) and ( p3| ) and d(u) in Eq. ( |B14| ). To obtain the 
series in e, we first expand the denominator in Eq. ( B28|) in powers of g*^ which is of order 
e. The computation is simple. The only integrals which require some manipulations are 
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+ — + — + - log 2- - bg 2 2 
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(2n) d v y m 2 + 4 
The final result is in agreement with our expression. 
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TABLES 



TABLE I. Three-dimensional estimates of g* from an unconstrained analysis, "unc", and con- 
strained analyses in various dimensions. For the analyses which use the estimates in d = 2 we 
report two errors: the first one gives the uncertainty of the resummation of the series, the second 
one expresses the change in the estimate when the two-dimensional result varies within one error 



bar. 


N 


unc 


d= 1 


d = 0,l 


d = 2 


d= 1,2 


d = 0,1,2 





1.37(9) 


1.39(3) 




1.392(23+1) 


1.390(16+1) 




1 


1.37(9) 


1.41(2) 


1.39(3) 


1.400(17+0) 


1.395(16+0) 


1.397(8+0) 


2 


1.36(7) 


1.39(2) 


1.41(2) 


1.401(15+2) 


1.411(12+3) 


1.413(8+5) 


3 


1.35(8) 


1.37(2) 


1.38(2) 


1.379(11+1) 


1.386(7+3) 


1.387(3+4) 


4 


1.33(5) 


1.35(2) 


1.36(2) 


1.357(11+3) 


1.364(11+6) 


1.366(7+8) 


8 


1.29(5) 


1.28(3) 


1.28(2) 


1.289(15+2) 


1.295(4+3) 


1.299(7+5) 


16 


1.20(2) 


1.198(13) 


1.196(10) 


1.198(8+0) 


1.198(4+1) 


1.199(1+2) 


24 


1.154(12) 


1.150(10) 


1.148(8) 


1.149(7+0) 


1.148(5+0) 


1.148(3+0) 


32 


1.123(10) 


1.121(8) 


1.119(7) 


1.120(6+0) 


1.118(4+0) 


1.118(3+0) 


48 


1.087(8) 


1.086(5) 


1.085(5) 


1.085(4+0) 


1.084(3+0) 


1.084(2+0) 



TABLE II. 


Two-dimensional estimates of g* obtained from analyses 


constrained at d = 1 and 


at d = 0, 1. 






N 


d = 1 


d = 0,1 





1.69(7) 




1 


1.79(5) 


1.75(5) 


2 


1.75(5) 


1.79(3) 


3 


1.68(6) 


1.72(2) 


4 


1.61(5) 


1.64(2) 


8 


1.45(5) 


1.45(2) 


16 


1.28(3) 


1.28(1) 


24 


1.21(2) 


1.20(1) 


32 


1.16(2) 


1.16(1) 


48 


1.11(1) 


1.11(1) 



32 



TABLE III. Estimates of /* obtained from the analysis of the available strong-coupling series 
with (GRT) and without (ST) the generalized Roskies transform (cf. Eq. (f74|)). We also report the 
values of /3 C and A used in our GRT analyses. The apparent uncertainty in the estimate obtained 
by employing the GRT is expressed as a sum of two numbers: the first number comes from the 
analysis at b = A, the second one is due to the uncertainty on A, and it is obtained by varying b. 
For N = 0, 1, 2, the estimate of /* is increasing with increasing b in the GRT. 



N 


lattice 




A 


GRT 


ST 





cubic 


0.213492(1) H| 


0.50(5) 


17.50(19+6) 


19.1(4) 


1 


cubic 


0.2216544(6) (7| 


0.50(5) 


23.64(14+10) 


24.9(3) 


1 


b.c.c. 


0.157373(2) @ 


0.50(5) 


23.53(9+10) 


24.9(2) 


1 


f.c.c. 


0.102062(5) [H 


0.50(5) 


23.53(18+10) 


24.6(4) 


1 


diamond 


0.36969(10) Q 


0.50(5) 


23.47(24+12) 


25.1(4) 


2 


cubic 


0.22710(1) 


0.52(2) 


28.45(17+5) 


29.4(3) 


3 


cubic 


0.231012(12) j7§ 


0.55(2) 


32.24(21+5) 


32.39(2) 


4 


cubic 


0.23398(2) || 


0.57(2) 


35.10(30+10) 


34.7(1) 


8 


cubic 


0.24084(3) || 


0.66(2) 


41.50(20+10) 


40.3(2) 


16 


cubic 


0.24587(6) || 


0.77(2) 


45.81(10+10) 


44.7(1) 


24 


cubic 


0.24795(3) 


0.83(2) 


47.28(6+10) 


46.4(1) 


32 


cubic 


0.24907(2) 


0.88(2) 


47.96(4+10) 


47.37(10) 


48 


cubic 


0.25023(2) 


0.93(1) 


48.66(4+4) 


48.32(10) 



TABLE IV. Summary of the three-dimensional estimates of g*. In Refs. 
were reported without errors. 



the results 



N 


e-exp. 


ff-exp. 


H.T. 


1/iV-exp. 





1.390(17) 


1.421(8) |||, 1.39 § 


1.393(20) 




1 


1.397(8) 


1.414(3) @, 1.40 § 


1.406(9) 




2 


1.413(13) 


1.406(4) |||, 1.40 [§ 


1.415(11) 




3 


1.387(7) 


1.391(4) §§, 1.39 [| 


1.411(12) 




4 


1.366(15) 


1.374 1@ 


1.396(16) 




8 


1.295(7) 


1.304 mnj 


1.321(10) 




16 


1.199(3) 


1.208 g|| 


1.215(5) 


1.204 


24 


1.148(3) 


1.154 ggg 


1.158(4) 


1.151 


32 


1.118(3) 


1.122 g|| 


1.122(3) 


1.1196 


48 


1.084(2) 




1.084(2) 


1.0839 
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TABLE V. Summary of the two-dimensional estimates of g* . 



N 


e-exp. 


H.T. 


1/N-exp. 


ff-exp. 


M.C. 





1.69(7) 


1.679(3) 








1 


1.75(5) 


1.7540(2) 




1.85(10) 101 


1.71(12) [281 

\ / 1 " 1 


2 


1.79(3) 


1.810(10) 






1.76(3) p| 


3 


1.72(2) 


1.724(9) 


1.758 


1.749(16) [Ml 


1.73(3) i9| 

i '~ J \ '~ J J II'-' ^|| 


4 

8 


1.64(2) 
1.45(2) 


1.655(16) 


1.698 
1.479 






16 


1.28(1) 




1.283 






24 


1.20(1) 




1.200 






32 


1.16(1) 




1.154 






48 


1.11(1) 




1.106 
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FIGURES 



FIG. 1. Plot of the function /(/?) of the three-dimensional Ising model as obtained from the 
quasi-diagonal PA's calculated with and without the use of the GRT (for each case we draw two 
lines representing the corresponding band of uncertainty). For comparison some Monte Carlo data 
taken from Refs. pH] (a) and [31] (b) are also plotted. 



FIG. 2. Two-dimensional 0(A r ) models with N = 3 and N = 4: Plot of /(/?) vs. £ as obtained 
by the analysis of its strong-coupling series. 



FIG. 3. Three-loop Feynman graphs contributing to the four-point correlation function. 
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